Edge spin accumulation in a ballistic regime 
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We consider a mesoscopic ballistic structure with Rashba spin-orbit splitting of the electron spec- 
trum. The ballistic region is attached to the leads with a voltage applied between them. We calculate 
the edge spin density which appears in the presence of a charge current through the structure due to 
the difference in populations of electrons coming from different leads. Combined effect of the bound- 
ary scattering and spin precession leads to oscillations of the edge polarization with the envelope 
function decaying as a power law of the distance from the boundary. The problem is solved with the 
use of scattering states. The simplicity of the method allows to gain an insight into the underlaying 
physics. We clarify the role of the unitarity of scattering for the problem of edge spin accumulation. 
In case of a straight boundary it leads to exact cancellation of all long-wave oscillations of the spin 
density. As a result, only the Friedel-like spin density oscillations with the momentum 2kF survive. 
However, this appears to be rather exceptional case. In general, the smooth spin oscillations with 
the spin precession length recover, as it happens, e.g., for the wiggly boundary. We demonstrate 
also, that there is no relation between the spin current in the bulk, which is zero in the considered 
case, and the edge spin accumulation. 

PACS numbers: 72.25.-b, 73.23.-b, 73.50.Bk 
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Currently there is a great interest, both experimental 
and theoretical, to spin currents and spin accumulation 
in various mesoscopic semiconductor structures [l|, |2]- 
Both phenomena are due to spin-orbit coupling and are 
of great importance for future spin electronics. The edge 
spin density accumulation related to the Mott asymme- 
try in scattering off impurities was recently measured [3| . 
Moreover, the edge spin density in the 2D hole system 
which is due to the so-called intrinsic mechanism of the 
spin-orbit interaction was also observed [4)] . One of the 
question of the interest is the relation between the spin 
current in the bulk and the spin density accumulated near 
the sample boundary. It is well known that in the diffu- 
sive regime (when a spin diffusion length is much larger 
than a mean free path) , such a relation indeed exists [3] , 
and the spin density appearing near the boundary is en- 
tirely determined by the spin flux coming from the bulk. 
For example, in the diffusive regime and in the case of 
Rashba Hamiltonian, when the spin current in the bulk 
is zero, the spin density component perpendicular to the 
lane is zero everywhere down to the sample boundary 
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However, in an opposite case, when the spin preces- 
sion length is much shorter than the mean free path, the 
situation is much less clear. An example of such a situa- 
tion is a mesoscopic ballistic (or quasi-ballistic) structure 
with spin-orbit-rclated splitting of the electron spectrum 
characterized by the energy A^, when one has the limit 
AflTp ^ 1, where Tp is the mean free time. Besides, in 
the presence of spin-orbit interaction, the boundary scat- 
tering itself is the source of generation of the spin density. 
It is quite clear that the characteristic length near bound- 



ary where the spin density appears is the spin precession 
length, Ls — Tivp/Aji, with vp being the Fermi veloc- 
ity. Thus in general there are two sources of the edge 
spin accumulation, the spin generation by the boundary 
scattering itself and the one due to the spin flux coming 
from the bulk (as it should be, at least, in the case of 
2D holes with a qubic s-o Hamiltonian). It is a sepa- 
rate problem to find the relative contributions of these 
two mechanisms to the edge spin accumulation. Now we 
study only the first mechanism when the spin density is 
generated upon the boundary scattering. 

For that we consider a mesoscopic structure described 
by Rashba Hamiltonian in the ballistic limit when a mean 
free path is much greater than the sample sizes. The bal- 
listic region is attached to the leads, and a voltage V ap- 
plied between the leads causes a charge current through 
the structure, as shown in Fig. [1] We stress that the 
electric field is absent inside an ideal ballistic conduc- 
tor, therefore the edge spin polarization appears not as 
a result of acceleration of electrons by an electric field 
and associated spin precession. Instead, the edge spin 
accumulation is due to the difference in populations of 
left-moving and right-moving electrons. Combined effect 
of the boundary scattering and spin precession leads to 
oscillations of the edge polarization. Thus, there is no re- 
lation between the spin current in the bulk, which is zero 
in the case considered, and the edge spin accumulation. 

We solve the problem by using scattering theory, with 
scattering states coming from different leads of the struc- 
ture and, therefore, having different occupations. The 
simplicity of the method allows to gain an insight into 
the underlaying physics. The problem of the spin den- 





FIG. 1: Left: Schematics of the boundary specular scattering 
in the presence of spin-orbit couphng. Plus and minus modes 
are shown for the same energy and the same wave vectors 
along the boundary. Right: Geometry of the system. Voltage 
V applied to the ideal leads causes a charge current through 
the ballistic region. 



oscillations acquire in addition smooth envelope function 
with the characteristic length Lg [g]. Therefore, strictly 
speaking, this phenomenon is different from a spin-orbit- 
related accumulation of the spin density upon boundary 
scattering. 

However, the case of a straight boundary appears to 
be a rather exceptional case, while in general the smooth 
spin oscillations with the spin precession length Lg re- 
cover, as it happens, for example, for the wiggly bound- 
ary, or for scattering off a circular impurity in a 2D elec- 
tron system [g] . This is a consequence of the fact that in 
higher dimensions the conditions of the unitarity of scat- 
tering take a different form, as explained below. In all 
these situations the spin density decays as a power law 
of the distance from the scatterer. 

The straight boundary. The s-o Hamiltonian in the 
bulk of a ballistic 2D electron system takes the following 
form 



sity accumulation in the ballistic case and for the straight 
boundary was erlier considered in Ref. [6| with the help 
of the Green's functions method (see also Ref. [7|). Sur- 
prisingly, in this case all the long- wave oscillations of the 
spin density cancel, and the final result contains only the 
Friedel-like oscillations with the momentum 2kp. Unfor- 
tunately, the method used in Ref. [6| does not allow a 
simple interpretation of the puzzling result. 

Here we show that it is the unitarity of scattering 
that leads to exact cancellation of the long-wave oscil- 
lations of the spin density with the Lg period in case of 
a straight boundary. It should be also mentioned that 
the observed behaviour is closely related to the effective 
one-dimensional character of the scattering problem aris- 
ing from the translational invariance along the boundary. 
This result may be interpreted as spin-orbit splitting of 
the Friedel oscillations in the charge density. In other 
words, two charge oscillations corresponding to spin- up 
and spin-down orientations, which coincided without s-o 
coupling, get shifted with respect to each other in the 
presence of the s-o interaction. Thus, short-wave length 
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where n is the normal to the plane, p is the 2D momen- 
tum, and M = ± are the helicity values. The solutions of 
this Hamiltonian have the form exp(ipr/?i)xAf (p), where 
r = x,y and the explicite form of the spinors is 
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with if being the angle between the momentum p and 
the positive direction of the x-axis. Now we consider the 
semi-infinite system and choose the x-axis to be directed 
perpendicular to the boundary (x = 0) of the 2D system, 
see Fig.[T] The wave functions which obey zero boundary 
conditions at x = arc obviously the scattering states 
which are the superposition of the bulk solutions and 
constitute the complete set of the orthonormal functions. 
Two scattering states corresponding to incident plus and 
minus modes with given wave vector along the boundary 
and the same energy are 
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The wave vectors are defined in the following way 
k^ = kl-kl kl=kl-kl, hk±^m{vFT^), (4) 



sin(y9 = ky/k^ and sin(/3i — ky/k^, as indicated in Fig.[T] 
From Eqs.([2]) and ([3]) one finds the scattering amplitudes 
F+ and F^: 



where hk± are the momenta at the Fermi energy in the 
plus and minus modes. The angles (p, ipi are defined via 



Fl 






Fl 



2cos(p 



(5) 



One can check that the amphtudes FZ and F^ for the 
incident minus mode with the same ky and the same en- 
ergy are obtained from i^^ and F^ by the replacement 
if -^r^ If I. The components of the unitary scattering ma- 
trix S have the following form: 



S^ — F^ , S_ — F_ , D_|_ — S_ — -F|_ 
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where Vx,i — dei{p)/dpx are the group velocities. 
Note that for the Rashba model one has Vx,-/vx,+ — 
cos (fi I / cos (fi. 

The wave functions Eqs. ([2]) and ([3]) may now be used 
to calculate the average z-component of the spin, Sz, as 
a function of coordinates: 



< Szix) >-- 
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where fpi^, ky) is the Fermi distribution function of the 
lead with the chemical potential shifted by ±eV/2. De- 
pending on the sign of ky , this function takes either of two 
values, /^(e — ^ — eF/2) or /i?(e — ^-l-eF/2). Substituting 
scattering states ([2]) and ^ to Eq. ([7]), we find that one 
may distinguish various contributions to < Sz {x) > with 
different oscillation periods, which originate from an in- 
terference of different terms in Eqs. ([2]) and ([3]). The 
part of < Sz (x) > which involves the interference of the 
outgoing waves (two last terms in Eqs. (I2]),(|3])), reads: 
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Here we used the fact that the distribution function 
fpi^, ky) has the same value at given energy for both plus 
and minus modes. Note, that the period of oscillations 
of the exponential factor e«(fe-fci)^ in Eq. ([S]) is of the 
order of the spin precession length. However, the term 
dH) vanishes, because the expression A is exactly equal 
to zero. Indeed, this expression is nothing but a non- 
diagonal component of the identity matrix SS'. Thus, 
we obtain an interesting result that the only reason for 
the absence of the long- wave oscillations in < Sz (x) > is 
the unitarity of scattering. 

Finally, by taking into account in Eq. ([7]) the terms 
responsible for the interference between incoming and 
the outgoing waves (for example, between first and sec- 
ond terms in Eq. [2]), and adding the contribution from 
the evanescent modes [10|, we reproduce Eq. (16) of 
the Ref. [6|. There it has been shown that so derived 
spin density may be approximated by the expression 




FIG. 2: Schematics of the scattering of the plus incident mode 
by the wiggly boundary, x = W sm{2TTy / X) . Besides the main 
scattering channels (solid lines), one has additional scattering 
waves with the wave vectors along the boundary shifted by 
±27r/A (dashed lines). 



Sz{x) oc eF cos(2fci?x) sin {x/Ls)/x. Therefore, the to- 
tal spin per unit length along the boundary scales as 
Jp dxSz{x) oc a^. Note, that the main contribution to 
this integral comes from small distances from the bound- 
ary, a; ~ Af, which makes it to depend strongly on the 
boundary quality. In what follows we show that, indeed, 
breaking the translation invariance along the boundary 
induces long wave length oscillations of the spin density 
and leads to the total spin that is not small in the pa- 
rameter a. 

Scattering by the wiggly boundary. We consider here 
scattering off the wiggly boundary, shown in Fig. [51 Be- 
cause in this case the translational invariance is broken, 
the condition of the unitarity of scattering takes a differ- 
ent form as compared to the case of the straight bound- 
ary, and, as a result, the cancellation of the smooth 
spin density oscillations does not take place. In order 
to demonstrate this effect we consider the mathemati- 
cally simple case of the abrupt impenetrable boundary 
described by the equation: x — C{y) = Vl^sin(27rj//A). 
To the lowest order in W the boundary condition reads: 



dx 
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We are looking for the solution in the perturbative form 
"^±{^,1/) = '^± ix,y) + '^j. {x,y), where the first or- 
der correction is proportional to W, and the zeroth 
order functions are given by Eqs.([2|),([3]). Using that 



*f^^(0,2/) = 0, we obtain 
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The following wave functions satisfy the above equation: 
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with the coefhcicnts a, 6, etc., being easily found from 
Eqs. dSJ, dSI, and ^. The definitions of the k- vectors 
are 



k^ =\lkl- {ky ± ^)2, k< ^\lkl~ {ky ± ^)2, 
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and angles are defined as (^ j^ ~ ipi + Sip^ ,(pf — ipi + Sipf , 
and so on. To the first order in l/\kp, 

Sipf = ±2n/Xki, 6f< = ±27r/Afc. (11) 



These equations are valid under the condition 
(27r/Afc_F) < cos^ ip. 

Calculation of < Sz >. To simplify somewhat long 
expressions, from now on we set 1/A = in the pre- 
exponential factors and keep only the terms of the lowest 
order in a. In particular, in view of Eq llll we disregard 
the difference between ip^ and ip, as well as between tpf 
and pi, etc., in the expressions for spinors. In other 
words, we assume Xp <S^ Ls <S^ A. Then, the first order 
correction for < Sz > in the minus mode reads 
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One can easily see that the contribution to < ^z > 
from the plus mode differs from the contribution from 
the minus mode by replacing everywhere A; o fci and 
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p -1^ pi. Thus we arrive at the following expression for 
the quantity in question: 
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Here the operation (A: ^ fci) should be understood as 
follows: fc — > fci, fcf — >■ fc>, etc. Note that this equa- 
tion takes into account only the propagating modes. It 
means that ky < fc+ , and for both plus and minus modes 
the wave vectors in x-direction are real. If, on the other 
hand, k+ < ky < k-, then for the given energy in the plus 
mode the wave vectors in x-direction are purely imagi- 
nary. Those modes are called evanescent. We will see 
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that in contrast to the case of straight boundary there is 
no exact cancellation of the contribution of the propagat- 
ing and evanescent modes into the edge spin density with 
^ oscillations (see below). We present here the expres- 
sion for < ^_|S'z|$_ > which is due to the evanescent 

modes in the case when all quantities k — \ k"^ ~ fci, 



K < = Jiky ± 2-k/\Y — fc+ are real 
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In order to obtain an analytical expression for the spin 
density, we expand Eqs. (TTSl) and (1141) with respect to 



2ttx/\, the actual small parameter being x/\/XXf'- 



J 



<s,ix,y) >=-C(y) 

(k + ki) 



dx 



.( ).W^cos(-^) 



{2Tr)^mvy A 




dkykll — cosffc — ki)x 

" y kki 



— cosffc + A:i)a; + 2cos2fca; + 2cos2fcia;l + / dfc„fc,^[2cos2fcia; — 2e^'*^(cosfcia; H sinfcia;)l>. (15) 

kki ^ ' Jk+ ^ K " j ^ ' 



r 



We note that Sstri^) describes the spin density distri- 
bution for the straight boundary, see Ref. [6|. There- 
fore, the first term in Eq. (jlSp reflects a small vari- 
ation of the boundary shape. The corrected solution 
Sstrix — C.{y)) contains only 2fci? -oscillations. Let us 
discuss now the next terms in EqllSI proportional to 
1/A. As it is clear from Eq. 1151 in general there are 
oscillations with three different periods, namely, 2kp- 
oscillations, and the oscillations with two long periods: 



^ and Ls- Here ^ = 1/i/fc?. — k^ is the new length 
scale, and under the conditions considered in the paper 
we have the set of inequalities Ai? <C ^ <C Ls, where 
Ls = l/(fc- — fc+) = fi/ma is the spin precession length. 
When ky —> fc+ (i.e. fc — ?► 0), then fci tends to 1/^, which 
clarifies the physical meaning of ^. 

We calculate now the integrals entering expression Eq. 

m 



h,2 - T- 



1 

kp 



ClrvyKy 



fcfci 



cos(fc =F ki)x, 



where the upper (minus) sign corresponds to /i and the 
lower (plus) sign to /2. Using new variables (see [uj), we 
obtain 



h = dz ■ \J z^ — [l/LsY cos{xz) = 

= ^^^^M) ^ £^^ + :^N,ix/L,), (16) 
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where Ni{x) is the Bessel function of the second kind, 
and 



h = 4fc|' / dzz v 1 — z^ cos{2kFXz). 
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Since l/2kp£^ ^ 1, we can separate the contributions of 
the long and short wave length oscillations for the prop- 
agating modes in Eq. [15] 
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Indeed, one can easily see that the integral l2kF con- 
tains only short wave length oscillations with 2kp mo- 
mentum, which may be disregarded. Therefore, the total 
contribution of the long wave length oscillations is given 
by 
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where the last term is the contribution of the evanes- 
cent modes. At the distances a:; ^ ^ we obtain the 
following dependence Itot = -2x^/3^^ -I- (a;/2L^)(7 -t- 
ln(a;/2Ls)). In the opposite limit, x 3> ^, we obtain 



Itot = ^Ni{x/L,) - i + 2£2£l£Zi)_ At even larger dis- 
tances, a; 3> is, one obtains smooth oscillations with 
the period of the order of Ls and the amplitude being 
proportional to ^/a. 

It should be noted that the term —1/x in Eq. (fT8|) 
has its origin in the interference of the incoming and 
outgoing waves, within each mode, with fc — >■ and 
fci — >■ 0. Those waves travel almost parallel to the 
boundary, and have the scattering angle close to zero. 
Note that the total spin per unit length along the 
boundary is proportional to the integral J dxltot {x) — 

{Tr/2Ls)J^dxNi{x/L,) - jf^ dx/x 
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is not small in 



-(l/2)ln(A/Af), i.e., it 
s-o coupling (which can be also seen directly from Eqs. 
(IT5)l and (HH)), in contrast to the case of the straight 
boundary. 

In conclusion, we have considered the problem of the 
edge spin accumulation in the ballistic mesoscopic struc- 
tures with spin-orbit-related splitting of the energy spec- 
trum. In the presence of the charge current the spin po- 
larization oscillates at the edge of the structure, as a re- 
sult of the spin precession upon scattering off the bound- 
ary. The scattering states method, being more trans- 
parent compared to the methods used earlier by various 
groups, has allowed us to clarify the physics associated 
with this effect. We have found that the boundary con- 
ditions play a crucial role, even the shape of the sample 
boundary being important. In particular, in the case of 
a straight boundary the absence of spin density oscilla- 
tions with the long-wave length, found in earlier theo- 
retical works, can be easily explained by the unitarity of 
scattering. On the other hand, the spin density accumu- 
lation in the case of a wiggly boundary has qualitatively 
different character as compared to the one for the case of 
a straight boundary, namely, the spin density oscillates 
with a large period of the order of the spin precession 
length, quite similar to the case of electron scattering 
off antidots in 2D system with the Rashba Hamiltonian 
[2]. Our analysis suggests that a smooth disorder at the 
boundary should induce spin oscillations with the long 
wave length. It also implies that in the ballistic case 
there is no relation between the spin current in the bulk 
and the spin accumulation at the edge of an electron sys- 
tem. 
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the absence of long-wave oscillations of the spin density 
coincide. This is no longer the case for 2D systems. In- 
deed, for a 2D system with the Rashba Hamiltonian the 
scattering problem of electrons off a circular scatterer 
(an antidot) may be solved exactly in the presence of the 
charge current. Then, the condition of the unitarity of 
scattering has completely different form as compared to 
the one for the case of scattering off the straight bound- 
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lating with period Ls, appears around the antidot and 
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When k+ < ky < k-, then for the reflected plus mode the 
wave vector in x-direction is purely imaginary. This mode 
is called evanescent. Of course, this situation only hap- 
pens for the incident minus mode, since only this mode 
can come from the leads under the conditions indicated. 
To calculate 7i we introduce the new variable z — 

ki ~ k = y fc^ + (1/5^) — k, which varies in the inter- 
val 1/Ls < z < 1/^. In terms of this variable the quan- 
tity ky = y/kl - fc2 = ^22Li - 1 + 2^^z2 ~ ^iz'^/2z£,^. 
Since ^ -C Ls, we always have \/k'^ — fc^ = 

k+y/z'^ - {1/Lsy/z. Note, that ^^fc^ = Ls/2. Finally, 
using that kydky — —kdk, and dz — —z-dk/ki, we obtain 
the equation of the main text. To calculate I2 we intro- 
duce the variable zi — ki + k — -v/fc^ + 1/^'^ + k, which 
varies in the interval 1/^ < 21 < (fc+ + fc_) = 2kF- Since 
^ <^ Ls, one gets for the quantity ky — \Jk^_i^ — fc^ — 

(l/2)-\/4fc|, — z\. Using that dz\ = z\ ■ dk/k\, and intro- 
ducing the variable 2 = z-i_/2kF, we obtain the expression 
given in the text. 



